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Abstract 
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■ We define the energy of a perfectly isolated system at a given retarded time 

as the suitable nuh hmit of the quasilocal energy E. The result coincides 
with the Bondi-Sachs mass. Our E is the lapse-unity shift-zero boundary 
value of the gravitational Hamiltonian appropriate for the partial system S 
contained within a finite topologically spherical boundary B = dT,. Moreover, 
we show that with an arbitrary lapse and zero shift the same null limit of the 
Hamiltonian defines a physically meaningful element in the space dual to 
supertranslations. This result is specialized to yield an expression for the full 
Bondi-Sachs four-momentum in terms of Hamiltonian values. 
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INTRODUCTION 



To define quasilocal energy in general relativity, one can begin with a suitable action 
functional for the time history Ai of a spatially bounded system S. Here "suitable" means 
that in the associated variational principle the induced metric on the time history T of 
the system boundary B = 9S is fixed. In particular, this means that the lapse of proper 
time between the boundaries of the initial and final states of the system E must be fixed as 
boundary data. The quasilocal energy (QLE) is then defined as minus the rate of change of 
the classical action (or Hamilton- Jacobi principal function) corresponding to a unit increase 
in proper time. [^,0 So defined, the QLE is a functional on the gravitational phase space of 
S, and is the value of the gravitational Hamiltonian corresponding to unit lapse function and 
zero shift vector on the system boundary B. Although other definitions of quasilocal energy 
have been proposed (see, for example, the references listed in [jl|), the QLE considered here 
has the key property, which we consider crucial, that it plays the role of internal energy in 
the thermodjTiamical description of coupled gravitational and matter fields. 

In this paper we define the energy of a perfectly isolated system at a given retarded 
time as the suitable limit of the quasilocal energy E for the partial system enclosed within 
a finite topologically spherical boundary.^ For our choice of asymptotic reference frame the 
energy that we compute equals what is usually called the Bondi-Sachs mass. |^,^] As we 
shall see, our asymptotic reference frame defines precisely that infinitesimal generator of the 
Bondi-Metzner-Sachs (bms) group corresponding to a pure time translation. We also 

show that in the same null limit the lapse-arbitrary, shift-zero Hamiltonian boundary value 
defines a physically meaningful element in the space dual to supertranslations. This dual 
space element, it turns out, coincides with the "supermomentum" discussed by Geroch. 
Our results are then specialized to an expression for the full Bondi-Sachs four-momentum in 
terms of Hamiltonian values. It is already known that when B is the two-sphere at spacelike 
infinity, the quasilocal and Arnowitt-Deser-Misner |10| notions of energy-momentum agree. 

Our results therefore indicate that the quasilocal formalism provides a unified Hamil- 
tonian framework for describing the standard notions of gravitational energy-momentum in 
asymptopia. 

Before turning to the technical details, let us first present a short overview of our ap- 
proach. Consider a spacetime A4 which is asymptotically flat at future null infinity and 
a system {w, R, 6, 0) of Bondi coordinates thereon. The retarded time w labels a one- 
parameter family of outgoing null hypersurfaces J\f{w). The coordinate i? is a luminosity 



^Hecht and Nester have also considered energy-momentum (and "spin") at null infinity (for a 
class of generally covariant theories including general relativity) via limits of quasilocal Hamiltonian 
values. Q Their treatment of energy- momentum is based on a differential- forms version of canonical 
gravity, often referred to as the "covariant canonical formalism." For pure bms translations our 
results are in accord with those found by Hecht and Nester, although at the level of general 
supertranslations they differ. We provide a careful analysis of the zero-energy reference term 
(necessary for the QLE to have a finite limit at null infinity), and this analysis is intimately 
connected with our results concerning general supertranslations. 



2 



parameter (areal radius) along the outgoing null-geodesic generators of the hypersurfaces 
N'lw). The Bondi coordinate system also defines a two-parameter family of topologically 
spherical two-surfaces B{w, R). It suits our purposes to consider only a single null hypersur- 
face of the family J\f{w), say N'lw^), the one determined by setting w equal to an arbitrary 
constant ly*. The collection B{w^,R) of two-surfaces foliates N{w^)^ and in the R ^ oo 
limit these two-surfaces converge on an infinite-radius round sphere 5(w*,oo). To stream- 
line the presentation, we refer to our generic null hypersurface simply as A/"; and we use 
the plain letter B to denote both the A/'-foliating collection B{w^:,R) and a single generic 
two-surface of this collection. Now, should we desire a more general A/'-foliating collection of 
two-surfaces, we could, of course, introduce a new radial coordinate R. For a fixed retarded 
time w = w^, the new two-surfaces would then arise as level surfaces of constant R. However, 
we shall not consider such a new radial coordinate, because the new two-surfaces would not 
necessarily converge towards a round sphere in the asymptotic limit. At any rate, we could 
handle such an additional kinematical freedom, were it present, by assuming that along each 
outgoing null ray R approached i? at a sufficiently fast rate in the asymptotic limit. 

Our first goal is to compute the QLE within a two-surface B in the limit as B approaches 
a spherical cut of X+ along the null surface A/", and to show that this result coincides with 
the Bondi-Sachs mass: 

M„Jw^ = lim / d'^xy/ae. (0.1) 

Here e = (k — k\"''^)/K is the quasilocal energy surface density with k = Svr (in geometrical 
units) and a is the determinant of the induced metric on B. Recall that k denotes the mean 
curvature of B as embedded in some spacelike spanning three-surface S. Since both B and S 
are embedded in the physical spa.ceiim.e Ai, we sometimes use the notation e\''^^ = k/n. Also 
recall that fcl'"' denotes the mean curvature of a surface which is isometric to B but which 
is embedded in a three-dimensional reference space different than S. Here we choose the 
reference space to be fiat Euclidean space E^, i. e. we assign a fiat three-slice of Minkowski 
spacetime the zero value of energy. |l|] Although a definition of the zero-energy reference in 



terms of fiat space is neither always essential nor appropriate |TT|, it is the correct choice 
for the analysis of this paper. 

In order to define k, we must select such a three-surface spanning B for each R value. 
(For a single B many different spanning three-surfaces will determine the same k. In fact, k 
is determined solely by B and a timelike unit vector field on B, which can be considered 
as the unit normal of a slice E. Thus, the continuation of S away from B is not needed; 
moreover, such a continuation of S might not be defined throughout the interior of Ai. 
Therefore, though we speak of choosing a S three-surface to span B for each R value, we 
are really fixing only a timelike unit normal vector field at B.) For generality, we leave the 
choice of spanning three-surface S essentially arbitrary at each R value, but we do enforce 
a definite choice asymptotically. Heuristically, as i? — > oo the S three-surface spanning B 
approaches an asymptotic three-surface which spans a round infinite-radius spherical cut 
of X"*" (see the figure). Our construction is, as expected, sensitive to the choice of asymptotic 
three-surface Soo- Said another way, the QLE depends on the fieet of observers at B whose 
four-velocities are orthogonal to the spanning three-surface at B. Therefore, one expects a 
priori the expression on the right-hand side of ( p.lj ) to depend on the choice of asymptotic 
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fleet associated with the two-sphere at X"*". The asymptotic fleet we choose corresponds to 
a pure BMS time translation: each member of the asymptotic fleet rides along d/dw. Note 
that, although d/dw is everywhere timelike in Ai (at least in the relevant exterior regions), 
the extension of d/dw to X"*" in a conformal completion j\4 of the physical spacetime A4 is in 
fact a null vector which lies in X"*". (While we occasionally flnd it clarifying to make reference 
to the concept of a conformal completion, we do not explicitly use conformal completions 
in this paper.) Therefore, heuristically, one should envision Sqo as a spacelike slice which 
becomes null asymptotically (see the flgure). 

This paper is organized as follows. In a preliminary section we write down the familiar 
Bondi-Sachs form p|,p!2[] of the spacetime metric as well as asymptotic expansions for the 



associated metric coefficients. We also introduce on Ai two future-pointing null vector fields 
k'^ and (do not confuse with the mean curvature k). Both vector fields point everywhere 
normal to our collection B{w, R) of two-surfaces, and /c^/'^ = —1. Next, we construct on A/" a 
timelike vector field := ^k^+l^ (equality restricted to Af), which in our analysis will define 
for each B along A/" a spacelike spanning three-surface S. In §11 we use the S three-surfaces 
determined by to define an unreferenced energy surface density el""''^ = k/n for each B 
slice of J\f and then examine the asymptotic limit oik/ n. In §111 we consider the asymptotic 
expression for the fiat-space reference density el""^ = k^"^ / but give the derivation of this 
expression in the Appendix. Next, we assemble the results of the previous two sections and 
prove the main claim (|0.1|) . In §IV we examine the "smeared energy surface density," which 
is the Hamiltonian value corresponding to an arbitrary supertranslation. We then specialize 
our result for the smeared energy surface density to express the full Bondi-Sachs four- 
momentum in terms of Hamiltonian values. In §V we examine the smeared energy surface 
density via the spin-coefficient formalism, and show that it equals the "supermomentum" 
of Geroch as written by Dray and Streubel. |[T^ The Appendix is devoted to a detailed 
analysis of the reference term. 

I. PRELIMINARIES 

In terms of a Bondi coordinate system the metric of our asymptotically fiat spacetime 
A4 takes the standard form [|7,12] 



g^.Ax^dx'' = -UVdw^ - 2UdwdR + aab{dx^ + W''dw){dx'' + W^dw) , (1.1) 

where a, b are B indices running over 6, 0. We assume the following expansions for the 
various metric coefficients above:0 



^Up to the A remainder terms, our expansions in the radial coordinate R coincide with those given 
by Sachs; however, we do not assume that the A remainder terms are necessarily expandable in 
powers of inverse R, an assumption which would be tantamount to what Sachs calls the "outgoing 
radiation condition." Recently, Chrusciel et al. have shown that "poly homogeneous" expansions 
in terms of log-'i? also provide a consistent framework for solving the characteristic initial value 



problem of the Bondi-Sachs type. |12| They argue that the so-called outgoing radiation condition 
is overly restrictive. 
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V = 


1 - 2mR-^ + Ay 


(1.2a) 


u = 




dUb) 


w' = 


{2X cot + X^g + CSC e)R-'^ + Aw9 




= 


CSC 6 {2Y cot + Ye- X <^ esc 6) i?"^ + A^^ 


(B) 


<^ab = 


RHab + \{2X)e^ae,b + (4r sin 9)e^(a<P,b) - (2X sin^ 0)</),,</), J R + A.„, . 


Od) 



Here X(w, 9, 0) and F(w, 6', 0) are respectively the real and imaginary parts of the asymp- 
totic shear c = X + lY , m(w, 9, (f)) is the all- important mass aspect, Sab is the metric of a 
unit-radius round sphere, and commas denote partial differentiation. In the Appendix we 
examine the form of the two metric aab in more detail. Remainder terms, denoted by the 
A symbol, always fall off faster (or have slower growth, as the case may be) than the terms 
which precede them. For instance. Ay denotes a term which falls off faster than 0{R~^). 

Introduce the future-directed null covector field fc^ = —e'^V^w, where the scalar function 
f] = rjiw, R, 9, 0) is a point-dependent boost parameter. The null covector fc^ is orthogonal 
to the spheres B{w, R), and the function rj gives us complete freedom in choosing the extent 
of fc^ at each point of any B two-surface. We shall find it necessary later to assume that r) 
falls off faster than 1/\/R on every outgoing ray. Also define another future-directed null 



vector field which is orthogonal to the B{w,R) and normalized so that k^l^ = — 1. As 
one-forms these null normals are 

fc^dx^ = -e'^dw (1.3a) 

l^dx^ = -e-'^UdR - ^e'^^UVdw , ([r|b) 

while as vector fields they are 

k'^d/dxf' = e^-^d/dR (1.4a) 

l^d/dx^ = e-^d/dw - ^e-'^Vd/dR - e-'^Wd/dx" . ([r|b) 



Now define := ^k^ + 1^^ and := ^k^ — l'^ along J\f as the timelike and spacelike unit 
normals of the B two-surfaces. For each slice B of the null hypersurface TV, the normals 

and n'^ determine a spanning spacelike three-surface E. As mentioned previously, the 
three-surface S is not unique and, moreover, need not be defined throughout Ai. Indeed, 
there is no guarantee that as defined is even surface-forming. (That is, in general does 
not satisfy the Frobinius condition U[Q,V^M,y] = 0.) Nevertheless, our construction provides 
us with what we need: a unit timelike vector orthogonal to B. We can therefore obtain 
an unreferenced energy surface density k/n which is the same for any slice or partial slice E 
that contains B and has timelike unit normal which agrees with sX B. 

Our construction implies 

u^djdx^ djdw (1.5) 

on each ray as i? — >■ cx). Now, the standard realization of the BMS-group Lie algebra (as 
vector fields on future null infinity) identifies the extension of djdw to X"*" (in a conformal 
completion of A^) with a pure time translation. Therefore, asymptotically, our fidu- 
cial surface Eqo determines precisely the pure time-translation generator of the BMS group. 
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We do not claim that u'^ generates an "infinitesimal asymptotic symmetry transformation" 
in the sense of Sachs 0, ^. e. that the various coefficients associated with the transformed 
metric (7^,^ + 2V(^'Ujy) satisfy the fall-off conditions ( |1.2| ); however, this is unimportant for 
our construction. 



II. COMPUTATION OF THE QUASILOCAL ENERGY SURFACE DENSITY 

We now turn to the task of calculating an expression for the unreferenced quasilocal 
energy surface density el^''^ = k/n. Our starting point is the definition k := —a^^V^n^, 
where the two-metric a'^'^ = g^" + 2k^^l'^^ serves as the projection operator into B. We find 
it convenient to write^ k = 2fi + p, where in the standard notation of the spin-coefficient 
formalism —fi and p are, respectively, the expansions associated with the inward null 
normal and outward null normal to B. These are given by the formulae 

/i = I^^^V^/a = |(V^P + k'^FV^l) (2.1a) 
P = -la^^V^k, = -i(V,.P + rk^v^k,) . d^b) 

As a technical tool, it proves convenient to introduce fiducial vector fields k^ and P deter- 
mined from ( |1.4| ) by setting = on M . From the middle expressions above, it is obvious 
that /i = e~'^jl and p = e^p, where the easier-to-calculate expressions fi and p are built exactly 
as in (|2.1|) but with the fiducial null normals k'^ and P. Therefore, we may assume that 1] = 



while calculating the spin coefficients in (|2.1|) and then simply multiply the r] = results by 
the appropriate factor to get the correct general expressions. Let us sketch the calculation. 
First, from ( p..3| a) with 7] = note that k^V^ku = 0, because k^, is a gradient. Next, using 
both expressions ( |L3|) with T] = 0, one can work the second term inside the parenthesis of 
(p?l|a) into the form k^l^V J.^ = —U~^l^V^U+]^Uk'^V^V. Finally, one writes the covariant- 
divergence terms as ordinary divergences; for example, V^A;^ = {—g)~^^'^d^{y/—gk^), where 
the square root of (minus) the determinant of the spacetime metric is y/—g = V 

Following these steps and multiplying by the appropriate boost factors at the end of the 
calculation, one finds 

p = ie-V-i(T - le-Wa-^a' - ^e-^SaW (2.2a) 
p = -ie^fJ-V-V. dJb) 



Here the over-dot denotes partial differentiation by d/dw, the prime denotes partial differ- 
entiation by d/dR, and 6a denotes the B covariant derivative. Since R is an areal radius, we 
may take a = R'^sin^ 6 [see the form of the B metric given in (|A4|) ]. Therefore, we obtain 
the compact expressions 

^ = -le-^VR-^ - le-^'daW (2.3a) 
p = -e^^u-^R-^ . d^^b) 



^Note that the definition of k does not depend on how is extended off B. 
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Adding twice (p73|a) to (pl3|b), we arrive at our desired expression 



(2.4) 



which has the asymptotic form 

k = -2R~^ + 2m(w„ e, 0)i?-2 - + Ak . (2.5) 

Note that we have chosen not to expand the 0{R~^) pure divergence term —SaW"". Our 
assumption about the fall-off of rj ensures that a term —1]"^/ R which appears in the asymptotic 
expression for k can be swept into A^. 



III. THE BONDI-SACHS MASS 

Write the total quasilocal energy as -E = — with the total unreferenced qnasilo- 
cal energy £'1^''^ taken as 

El^*^^ = - / d^xy/^k. (3.1) 

K Jb(w*,R) 

Plugging the expansion ( p.5|) into the above expression, using ^ = R^ sin 9 for our choice 
of coordinates, and integrating term-by-term, one finds 

E]^"^ = -R + M^^{w,) + Aei^^y. (3.2) 

Here the Bondi-Sachs mass associated with the w = w^, cut of X"*" is the two-surface average 
of the mass aspect evaluated at w = w^,, 

M^^.{w,) = - fdQm{w,,e,(j)). (3.3) 



K 

We use the notation / dfi := Jq d6 Jq^ d(j)sm6 to denote proper integration over the unit 
sphere (which is identified with a spherical cut of X+). In passing from ( |3.1| ) to ( p.2|) , we have 
made an appeal to Stokes' theorem to show that the "dangerous" 0{R^) term that arises 
from proper integration over the pure-divergence term SaW"" in ( p.5| ) does indeed vanish. 
Hence, this term does not contribute to the Bondi-Sachs mass and does not spoil the result 

The reference point contribution to the energy is 

= -- [ d'^xV^kl'^^' , (3.4) 

K Jb{w,,R) 

where the asymptotic expression for k^"^ must be determined from the specific asymptotic 
form ( [X^ ) of the Sachs two-metric. We present this calculation in the Appendix. The result 
is 

~ E^"' = R + Q ■ R^ + AE\r.f . (3.5) 
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Note the absence of an 0{R^) term in El""^. The result (3.5) has just the right form, in that 
it removes the part of i?!""'"^ which becomes singular as i? — oo but does not itself contribute 
to the mass. Therefore the total quasilocal energy for large R is 

E = f d'^x^/^e = M^^{w,) + AE. (3.6) 

Jb{w,,R) 

This is the energy of the gravitational and matter fields associated with the spacelike three- 
surface S which spans a B slice of A/" and which tends toward Sqo. Our main claim ( p.l|) 
follows immediately from (3.6). 



IV. SMEARED ENERGY SURFACE DENSITY 

Consider the expression Hb for the on-shell value of the gravitational Hamiltonian ap- 
propriate for a spatially bounded three- manifold S, subject to the choice of a vanishing shift 
vector at the boundary 9S = B: 

Hb= I d^x^Ne. (4.1) 
Jb 

We refer to Hb as the smeared energy surface density. Addition of this boundary term to 
the smeared Hamiltonian constraint ensures that as a whole the sum is functionally differen- 
tiable. In this section we consider the R oo limit of Hb along the null hypersurface Af 
in exactly the same fashion as we considered the limit (1.1) of the quasilocal energy previ- 
ously. Before evaluating lim/j_^oo Hb(w*,r), let us discuss its physical significance. Consider a 
particular spherical cut B{w^, oo) of . A general BMS supertranslation pushes B{w^,, oo) 
forward in retarded time w in a general angle-dependent fashion. As is well-known, the 
infinitesimal generator corresponding to such a supertranslation has the form ad/dw\j+, 
where a{9, (p) is any twice differentiable function of the angular coordinates. As we have 
seen, d/dw is heuristically the hypersurface normal at B{w^,, oo) of an asymptotic span- 
ning three-surface Sqo. In other words, each member of the fleet of observers at B{w^, oo) 
rides along d/dw. Therefore, again heuristically, the on-shell value of the Hamiltonian gen- 
erator of a general BMS supertranslation is 

/ d'^xy/^ae . (4.2) 

Jb{w,,oo) 

This symbolic expression coincides with the R oo limit of the smeared energy surface 
density (4.1), where we set a{6, (p) := limij_>oo N{R, 6, (p) (suitable fall-off behavior for is 
assumed). Thus, lim^j^oo Hb(w,,r) defines a physically meaningful element in the dual space 
of general supertranslations. In this respect it is like the "supermomentum" of Geroch. |^ 
In the next section we show explicitly that, in fact, lim^^oo Hb(w,,r) is precisely Geroch's 
"supermomentum." Note, however, that it might be better to call such an expression the 
"superenergy," as it arises entirely from the "energy sector" of the Hamiltonian's boundary 
term (that is, the sector with vanishing shift vector) but also incorporates the "many- 
fingered" nature of time (that is, an arbitrary lapse function). 

Let us now evaluate the R ^ oo limit of the smeared energy surface density Hb- As we 
have stated, there is no 0{R^) contribution to The absence of this contribution stems 
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from the fact that the two-sphere average of the coefficient ^'^''k\'"'^ of the 0(R~'^) piece of the 
reference term vanishes. As spelled out in the Appendix, this fact follows directly from 
an equation governing the required isometric embedding of B into Euclidean three-space. 
Moreover, as seen in Section 2, the coefficient ^"^^k of the 0{R~^) piece of the physical k is 
not solely twice the mass aspect but also contains a unit-sphere divergence term. Now, in 
the present case e = (A; — A;|''°')/fi; is smeared against a function A^, so one might worry that 
the limit is spoiled in some way by the presence of the smearing function. However, as we 
now show, for solutions of the field equations, the unintegrated expression AnR'^e is precisely 
the mass aspect of the system in the R ^ oo limit. This striking result rests on an exact 
cancellation between ^"^^k]""^ and the aforementioned unit-sphere divergence part of ^'^^k. 

With the machinery set up in the previous sections and the Appendix [see in particular 
equations (p.5|) and (A2)], we find the following limit: 

hni^ ^kR'^e = m(w„ ^, 0) - ^ [esc 6 da{sm O'^^W^) - ^'^^>r\ . (4.3) 

Here we set k = 8tc (in geometrical units) and use '•^•'T^ to denote the coefficient of the 
0{R~^) piece of the B Ricci scalar. Also, the coefficients (^W" of the leading 0{R~^) pieces 
of W"' are listed in (1.2c). Inspection of (|1.2| c,d) shows that ^'^^W'^ is expressed in terms of 
the same functions, X and Y, that appear in the 0(-R~^) piece of aab/R"^- Furthermore, a 
short calculation with the B metric shows that for these solutions may be expressed in 
terms of ^"^^W"" as follows: 

-i(3)7^ = _csc^(9„(sin^(2)y^''). (4.4) 
Therefore, the term in (|4.3[) which is enclosed by square brackets vanishes, and we obtain 



lim / d^x^Ne = - f dQa{e,(j))m{w,,e,(j)) , (4.5) 

R-^oo Jb(w,,R) hi J 

for the desired limit. This result shows that the R ^ oo limit of the smeared energy surface 
density equals the smeared mass aspect. Coupled with the findings of the next section, it 
follows that Geroch's "supermomentum" is just the smeared mass aspect. This simple result 
does not appear to be widely known. 

Finally, recall that the Bondi-Sachs four- momentum componentsQPj^ correspond asymp- 
totically to a pure translation. In terms of the smeared energy surface density, one ob- 
tains a pure translation for a judicious choice of lapse function on B{w^, oo); namely, 
a{9, (p) = e\a-{9, (p), where the ex are constants and ^ 

= 1 (4.6a) 

a- = sin 9 cos (j) (|4.6| b) 

a- = sin 9 sin ( [4.6| c) 

a^ = cos^. §3d) 



Therefore, we write ^xP^g{w*) = ^iT^R-*oo Hb{w,,r) for the appropriate limiting value of A^, 
and thereby obtain the Bondi-Sachs four-momentum as a Hamiltonian value. 



^Underlined Greek indices refer to components of the total Bondi-Sachs four-momentum. 
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V. SUPERMOMENTUM 



In this section we show that the null limit of the smeared Hamiltonian boundary value, 
Eq. (4.1), is the "supermomentum" of Geroch. ^ To be precise, we show that in the null 



limit Hb equals Geroch's "supermomentum" as written by Dray and Streubel. [13| The 



spin-coefficient formalism is required for this analysis.^ Apart from a few minor notational 



changes we adopt the conventions of Dougan. |]14| Geometrically, the scenario is nearly the 
same as the one described in the previous sections. However, we now work with a slightly 
different type of Bondi coordinates. Namely, {w,r,(,0, where r is an affine parameter 
along the null-geodesic generators of Af and ( = e^'^ cot(^/2) is the stereographic coordinate. 
Dougan pick^ = — V^w as the first leg of a null tetrad, which is the same normal as 
given in ( |1.3| a) if rj = 0. For convenience, in this section we ignore the kinematical freedom 
associated with the rj parameter, setting it to zero throughout. As before, the vector field 
^k^ + (equality restricted to J\f) defines a three-surface S spanning each B slice of 
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A/". It follows that u^d/dx^ d/dw as r ^ oo, and hence our asymptotic slice Sqo again 
defines a pure BMS time translation. 

Let us first collect the essential background results from which we will need. First, 



the required spin coefficients have the following asymptotic expansions:[] 

p = -r'^ - a^a\-^ + ©(r^^) 

a = a^'r-^ + 0{r-^) , (5.1) 

where is Dougan's — p', the term \E'2 is a certain asymptotic component of the Weyl tensor, 
and cr° is the asymptotic piece of the shear. Like before, an over-dot denotes differentiation 



by d/dw. As fully described in |14], is the standard differential operator from the com- 
pacted spin-coefficient formalism, here defined on the unit sphere. The expansion for the 
corresponding full operator in spacetime is 



^Throughout §V we deal exclusively with smooth expansions in inverse powers of an affine radius, 
as we know of no work examining the standard spin coefficient approach to null infinity within a 



more general framework such as the poly homogeneous one. The expansions we borrow from |14] 
are valid for Einstein-Maxwell theory. 

^Our and respectively correspond to la = V^n and in |14|, where u is Dougan's retarded 
time. The minus sign difference between our definition for k^ and Dougan's definition for la stems 
from a difference in metric-signature conventions [ours is (—,+,+,+)]. The convention for metric 
signature does not affect the spin coefficients (|5.lD. 



^Note the dual use of a as both the stem letter for the B two-metric and as the spin coefficient 
known as the shear. We have used a twice in order to stick with the conventions of our references 
as much as possible. In all but equation ( ^.ID , where a has the spin-coefficient meaning, it carries 
a "0" superscript denoting the asymptotic piece. 
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^ = + r-^[s(^oa°) - <y%] + 0(r-3) , (5.2) 

where s = sw{ip), ip being the spacetime scalar on which ^ acts and SW denoting spin weight. 
The commutator of ^ and ^ is 

(^^ - #)(^ = is7^(^ . (5.3) 

Now consider the following ansatz for the B intrinsic Ricci scalar: 

7^ = 2r~^ + (3>;^r^^ + 0(r-^) . (5.4) 

If we insert this expansion into ( [5.3| ) and expand both sides of the equation [assuming 
ip = i^)(p + (^^(^r~^ + 0(r~^) with sw{ip) = 1], then to lowest order, namely 0(r~^), we get a 
trivial equality. However, equality at the next order demands that 

i^^yjz = -^y + . (5.5) 

This will prove to be a very important result for our purposes. Finally, Dougan gives the 
following expansion for the B volume element: 

d^x^ = dnr\l - aVr-^) + 0{r-^) . (5.6) 

(Here a is the determinant of the B metric and is the asymptotic piece of the shear.) 

We now consider the spin-coefficient expression for the smeared energy surface density 
introduced in §IV. Again with k = 2fi + p, in the present notation we find 

k = -2r-^ - 2[^° + a^cT^ + f.cry-^ + 0{r-^) . (5.7) 

Moreover, by an argument identical to the one found in the last paragraph of the Appendix 
(although here with the affine radius r rather than the areal radius i?), we know that the 
result (|5.5|) determines 

= -2r-i - (^2^° + ^yy-^ + 0(r-3) (5.8) 

as the appropriate asymptotic expansion for the reference term. Therefore, (k times) the 
full quasilocal energy surface density is 

KB = -2[^° + + - l^^cry-^ + 0(r-3) . (5.9) 

At this point we consider again a smearing function A^, with appropriate fall-off behavior 
and limit a(C,C) = ^^^r^oo N{rX,C)- Using the results amassed up to now, one computes 
that the limit of the smeared energy surface density is 



\im f d'x^Ne = -l[dQa [^^ + a'cr' + \fy - \fy]\ . (5.10) 



2 

d^x^yaNe = — 

The right-hand side of this equation is the "supermomentum" of Geroch as written by 
Dray and Streubel [see equation (Al.12) of [|13| and set their 6 = for a Bondi frame as 
we have here]; and the "supermomentum" is known to be the "charge integral" associated 
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with the Ashtekar-Streubel flux of gravitational radiation at X+ (in the restricted case 
when the flux is associated with a supertranslation). |16| Dray and Streubel have discussed 
the importance of the particular factors of | which multiply the last two terms within the 
square brackets on the right-hand side of equation ( ^.101) . It is evident from our approach 
that the origin of these | factors stems from the flat-space reference of the quasilocal energy 
(flat-space being the correct reference in the present context). When a determines a pure 
BMS translation, the last two terms in the integrand integrate to zero. For instance, setting 
a = 1, one finds that the strict energy 

E = I d^x^e = -- [dn + (5.11) 
is the standard spin-coefficient expression for the Bondi-Sachs mass Mgg{w^). |T^JT7[] 
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APPENDIX: SUBTRACTION TERM 

In this Appendix we prove that the subtraction-term contribution —E\"'^ to the quasilocal 
energy obeys 

- E[^' = R + 0-R^ + A^|rof (Al) 
for large R. Moreover, we derive the result 

(2)A;|-f = _i{3);^ (A2) 

relating the coefficient ^'^^k\"°^ of the 0{R~^) piece of the reference term k]""^ and the coefficient 
of the 0{R~^) piece of the B Ricci scalar TZ. 
The —El""^ term is constructed as follows. Consider a generic B slice of A/". Assume that 
B may be embedded isometrically in Euclidean three-space E^ and that the embedding is 
suitably unique (we address these issues below). Let {k\'"'^)ab represent the extrinsic curvature 
of B as isometrically embedded in E^. The flat-space reference density is e\'"'^ = k\"'^ / and 
in terms of this density 

-E|-f = _/ d?x^er . (A3) 

JB{W:,,R) 

Since it is the intrinsic geometry of B which determines the reference term k\''°\ let us 
first collect a few results concerning this geometry. In the Bondi coordinate system the 
two-metric of B takes the following form: ^ 
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aabdx"-dx^ = 

\R^{e^^ + e^^)dd^ + 2i?2 gin 6 sinh(7 - 5)ded(t) + \R^{e-^^ + e'^^) sin^ ^d^^ . (A4) 
In terms of the functions X and Y , 7 and 5 have the expansions [see ( |1.2| d)] 

7 = (X + + 

(5 = (X - + A5 . 



(A5) 



As mentioned, we do not necessarily enforce Sachs' "outgoing radiation condition" which 
would, in fact, imply that in the expansions for 7 and 5 the terms that appear after the 
leading 0{R-^) terms are 0{R-^). [§ 

From the form of the line element ( |A^ ) one easily verifies (|1.2| d) and that ^/a = R^ sinO. 
A bit more work establishes that the scalar curvature TZ oi B has the asymptotic form 

7^ = 2R-^ + (^>/^i^-^ + At^ , (A6) 

where ^^-"T^ = ^^^{w^^O^cj)). We have given the explicit form (|4.4]) of the coefficient '•^•*7?. 



corresponding to the asymptotic solutions considered here, and we have found this coefficient 
to be a pure divergence on the unit sphere. That this term integrates to zero can also be 
shown via the argument to follow, which does not assume that the Einstein equations hold. 
Lemma: The two-sphere average of (3>7^ vanishes, i. e. J dQ ^^>}Z = 0. To prove the lemma, 
start with the Gauss-Bonnet theorem |18[ 



57r = / d^x^/aU 

= d^x sin ^ (2 + ^^>RR'^ + R^An) ■ (A7) 

In the second line we have simply expanded TZ and used y/a = R^ sin 6. On the right side 
of the equation, integration over the first term inside the brackets gives Sir. Therefore, we 
arrive at 

= J dQ^^>R + R^ J dQAn. (A8) 

Since / dQ A-ji falls off faster than 0{R~^), the R ^ 00 limit of this last equation proves 
the lemma. 

Now let us discern the R dependence of the reference term kl""^. We must first address the 
issue of whether or not B may be isometrically embedded in in a suitably unique way. 
It is known that a Riemannian manifold possessing two-sphere topology and everywhere 
positive scalar curvature can be globally immersed in E"^. An immersion differs from an 
embedding by allowing for self-intersection of the surface (seemingly allowable provided that 
^|rcf 2;g]2iains well defined). The Cohn-Vossen theorem states that any compact two-surface 
contained in E^ whose curvature is everywhere positive is unwarpable. Unwarpable means 
that the surface is uniquely determined by its two-metric, up to translations or rotations in 
E^ (which of course do not affect A; I""'). |1|,|T^ Our surface B has two-sphere topology, and 
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for R sufficiently large the scalar curvature TZ is everywhere positive. It thus follows that 
for suitably large R, the reference term k\'"'^ is well defined. 

The key equation for determining the form of /cl""' is the standard Gauss-Codazzi relation: 



{kn'-imkn':-T^ = o- m 

This equation (essentially a two-dimensional version of the Hamiltonian constraint) is an 
integrability criterion, obeyed by our embedding, which relates certain components of the 
vanishing Riemann tensor of to the intrinsic B curvature scalar and the desired reference 
extrinsic curvature tensor {k\"''^)ab- Next, since B approaches a perfectly round sphere as 
i? — > cxD, we take the following expansions for the various pieces of (^1""')^ as an ansatz: 

{kne = + ^MnlR-^ + A(fc) (AlOa) 

{knt = -R-' + ''^{kntR-' + A(.) imp) 
(A;|-o^ = o-i?-^ + A(,). mm 

To avoid clutter, for the various remainder terms we have used in ( |A1(]| ) simply a subscript 
ik) in place of what should be (fcl'^Oe^ ^^c. Plugging the expansions ( [A6| ) and ( |A1CI| ) into 
(|A9|), one finds that to lowest order, namely 0{R~^), the equation is identically satisfied. 
At the next order, namely 0{R~^), the equation ( |A9| ) yields the result (A2). Therefore, we 
have found the following asymptotic expansion for k\"'^\ 

^|.of _ _2^-l _ li^yjlR-^ + Afc|.ef . (All) 

The first lemma we proved above has an important consequence. It ensures that the "dan- 
gerous" 0{R^) term in the integral (|A3|) vanishes (regardless of whether or not ( [4.4|) holds. 



which requires that the Einstein equations are satisfied asymptotically). Hence we get the 
result (|AT1) . 
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FIG. 1. Geometry of the asymptotic slice Soo 
In this figure one dimension of the two-surfaces B{w^, R) is suppressed. The shaded, partially 
cut-away, conical surface depicts the null hypersurface Af = ^/{w^t ) determined by a constant value 
w = w* oi retarded time. Heuristically, in the limit i? — >■ oo the S slice spanning B{w^,R) becomes 
the asymptotic slice Soo which spans a round spherical cut of T"*", and one should envision the 
spacelike slice Eqo as becoming null asymptotically. Although d/dw is timclikc everywhere in the 
physical spacetime M (at least in relevant exterior regions), the extension to X+ of d/dw (in a 
conformal completion Ai of Ai) is a null vector which lies in Again heuristically, on B(wi,, oo) 
the S(X) hypersurface normal u^d/dx^ is d/dw. 
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